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PREFACE 


This document reports the preliminary processing and analysis 
efforts on one phase of a project designed to improve the classi- 
fication of forest inventories. The work focused on the differen- 
tiation of forest type -jignatures in hilly and mountainous 
terrain where the effects of low Sun angle cnduced a wide 
variation in like signature response. 

The authors wish to acknowledge the assistance of Dr. B. P. 

Edwards of Lockheed Electronics Company, Inc. for his editorial 
comments and recommendations pertaining to the illumination 
model. 

This type III document has been approved by the Supervisor of 
the Forestry Applications Section for distribution to persons 
directly assocxated with the Nationwide Forestry Applications 
Program. 
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1. INTKOlHk'TlON 

This i)ap«‘r ili*su*r iJ)os ihu procoduros and rt'sults of a vary bri(*t 
study ol the rolationsnip bidwaan iorast ralUu’tance, torrain 
s)o})i.‘ and a lov» ilun anaLa and its idioct on l»indsat spoatral data. 

To t‘xptdiU! the study, s impl i t lewl teclm viiut'S and approximat ions 
won' t'mployed and which ar(‘ deseriba*l tully in aj>pand*lx h. The 
study was stimulated by previous work conducted dui'iaif the Ten- 
Kcosystem ,‘5?udy (Tlk!) tor Washin-jton County, Missouri. Durinu 
t lu‘ photographic anaJysis phasi' ol tdiis pro'p'cl , it was observed 
that both the hidh altitude photoitraphy and bandsat imagery dis- 
pliiyed an unusual amount; ol shadow caus(?d by the hilly terrain. 

This eltect. was a.ccentuat «'d l>y the intensity of c,he f>un*s irradi- 
ance on the sLopo* lacing the soutlu>ast. (tlinnit ion of the Sun's 
TMvs) and the opiiosite slope wlu‘ri> the irradiance just yravied thi' 
surlaci*. Th*! rclativfi in* ight.aess ol tlu' two surfac’os gave rise 
to the initial assumption of shadows lieiiut projected by the 
terriun. This phenomenon was especially not.iceablci in the north- 
wi'Stern portion ol the county v.h<*re tlu^ tie rain was well disnetod 
by numoi ous headward working st roviras givimt a topographic expres- 
sion in the term ot narrow ridifos and poorly defined valleys. A 
later analysis ol the topographic data showed that the average 
slope was not stoop enough to cast a true shadow and that the 
observed plumomenon was caused by a combination of slope, Hun 
elevation, and iiidiroctional reflectance of the surface features 
(fig. 1). 

2. OBJKCTIVIH 

The objectives of this study are, two-fold); 

a. To determine the ef'fofts of the low Hun angle and slope on 

spectral/onvironmontal signatures obtarnod from Landsat imaijory 

* 

A shadow cast by an opaque body as compared to an apparent shadow 
(psGudoshadow) which may be attributed to the bidirectional charac- 
teristics of the surface or angle of illumination. 
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Figure 1.— Pseudoshadow effects in hilly terrain. 


b. 


Tv) t uul .1 invans ot ivdui'ijui thv' v'l tovle of Uii!'. phenomenon 
HO a« to im’rvMSt' the av’cnraf’y of lutuie computer 
cl4B«i licat ions 


i. HTUDV AKL'A 

Tin* Btiuiy ai vm, »h fentcr»Hl in t lu» nurt hwestt*rn vp'.adrant of 
Washnuiion (’onnty, MiHHtmii (tiq. 2). The toi)Ov|raphy of the area 
can la'st ho ilesc! ibi>d .ui consiBtiiuj ol moderately steep ridges 
with a mean uiavlieid u1 ajjproximiitely 20 d agrees and a maximum 
ijia»liv'nt ol U) vlvsptn'B, Tlu' axial trend ol tlu' slopes is 
iu>rt hwv':;t /'novtiluM:.! ; however, thin in Hub 'iin't to local variations. 
Thv* ri'iivd in appr»»ximat ely 77.2 meters (250 fv'v't) while the 
highv'Mt v‘l»'vatlon in l In* area in I’l int Hill, witli an elevation of 
402.1 meter;; (1120 leetj above mi*an stt'a li'Vv'i. 

The v(*q**tat ion <*o.. *st.*i of mixv'vl luirdwoods principally oak and 
oak”hii'kory , with minor stands of pim* and (’■edar. The pine and 
eediir plantativ)ns ocevr in Isolateil sfands and are usually 
v’onfini'il to tlu* u})pv‘r slopo!;. Due to the hilly tv*rrain and 
thiek foresf (’anopy, grans lands ar(* invarialdy locuaiod in the 
narrow val h*yfi and, for Uu* most i^ari, vronsist of well developed 
pantures. In fhis quadrant of the v’ounty, .ipproxiraately 85 per- 
cent of flu* art'a ir; li.irdwood, 10 perc(*nt grassland, and 5 percent 
softwood. 


4. MATKRIAh AND KQUIDMKNT 

The General Klv»ctric Imatje 100 lnierac?tivo Computer was utilized 

extensively for the image ilisplay vrapabilities of the system, 

signature aciiuisition, and classification. Landsat data from 

1 

scene 1845-1559 dated November 15, 1974' , was vitilized during the 

^This Landsat scene was utilized during the TKS when the 
phenomenon was first observed. 
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Figure 2. 


- Solo<:torl stwly sites In Wasliinqton County, Missouri 

(^^calo a})proximatoly 1 : 25 ( 1 , 500 ) . 





St miy a.s iiisj<layu.i *i t i-i h I'h.ti ti t t-t i :,t n* of a low Sun 

ari'il'* {JM litMji j . Hiqli itii.itaii- ^ol*>i iutr.it ml (iMK) i4ioto- 
«ji .ilihy : mm .'H*i , t iuwii ‘i* .* • sti. •• IT74, w.is utilizml as an 

aivi lu iuj 1 h' st-vi:. • i iu:.t at i.i. ta itiiMtiuu of v«.‘i}itt at i«m, 

aiui tat t li' it'.. I .-iit Miu at ihi< • -toliv with aoutt)!!!' data from 

tha 5*. f ‘,1 I .|;a.il Stu'iy }>' i ai*w< ml;; t.iua.lr.inq It* at l;<di OIH) 

sc'.ilu. !u oi d« I !i< vai ity t h*^ |i ‘i,.i I i ou aov«‘r ohscurod by t,ht* 

psi'udowsiiaaow.i, *!-•• toil(,w!a>} a J ati-j; durnvf a pel iuil of 

hi‘iht*r ,*:uu w.i:; it.ili. a i; h.ui i.'i.it tia.-ua 17 July fO, 

17 / 4 , and hi*;h altiiuda an* phof mm* tpuy *im;. “i.stnon t l«>v?n 

Juj>tuinh(‘i l't/4. 
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The.* .ippro.ich waa :J ! .i i lilit t at wat;»i - spat-i*,*! u!* atjuramnnts of tho 
Laadsat hiiaMa.'y wata oldaiu',‘1 of 'ialaat.ii .na,*;: whU’h disixlayed 
well tied lauil .ajaAtratit slMdawi utj aHiMdo, .iiii .similar, if not 
identii’al, voaatai.iou aovat .tuu as];ai*f .tuela. t’oireetion factors 
wen* then ' biaiuad lot is.n’h sLope* .md .ipjd lad to the initial 
moasuremr*nf !j, thus <*oito>*lJuj .d 1 lata t-i .i !'.uu .inqlct normal to 
the surla'T*. M wan hypoth* *: j i ;tad that th<‘ a|>pi Laaiion of tho.se 
corr ections to the h.adwood s i jnaf.urt^s would '’itiluee the variation 
in spectral resjornse ij tha satnatures i^ilJained from the slopes 
for all .mylo.s. Assuming thi.*} cjouI'I 1.h,* ar*hiavr*(l, with only 
minor residaal.s, Uu-n means gouM l?e davalopi-d whereby these 
corrected vaLue.s could be ai>plied to Uu* Jiaiidsat’. data with the 
intent of improviiij t lu* iiai.i prioj to r* i.e;*; j f ieation. It should 
be noted that; the el toid.'; of al mospluu ia f .iators, chanyc of 


aspect, bun a/.imui.h, sensor look anjb* (and a host of others) 
have not Ijoeu eons i'ltired. as jiarf of thi.*; reporf since f he main 
theme is direcfaMl f.owards inaint ainiuj a c.inipl l.stic apjiroach. 


b. 1‘PKMMlNAKy TLCIlNlUw. . 

For thu purposo of ^HueBsinq iht* of feints of tho apparent tihad,ow* 
or "i soudoHhaciowH" an shown in f iqiiro 1, sovon potential stud' 
situs in the foriti of transects were selected fiom the Landsat 
ima<jety« Kach transect was characterized oy certain features 
which niiqiit influenct* the spectral lesiions**, i.e., azimuthal 
orientation of the transect, vegetation cover types, slope and 
aspect, and the }iscudoshadow id feci s. Care was exercised to 
insure that the veifi*tation cover alomj tlu? transect was uniform 
and that each transect was oric*nted aormal to the mean axial 
trend of the ridqea so that the Landsat radiometric measurements 
could be recorded along the diiection oi maximum slope. Prior 
to mensuration of the rad-iince values, the transects were trans- 
ferred to the aerial phoiugraphy so that the terrain slopes 
could be accurately identified and correlated with the radiance 
values. Once the slopes woie idwitiJ led, the aspect angle was 
then determined by transposing the trace of Lhe transect to the 
USGS map source ,ind computing the detjrec ol slope from the 
contour information and base itistance of the slope. 

To implement the .study, a .spectral re.soon.se profile (fiff. 3) was 
made across Tran.sect B (fig. 4) to provide preliminary data to 
determine how the s}»ectral response varied with the topography 
and vegetation cover. The spectral measurements were obtained 
from the Landsat imagery utilizing the Image 100 system. The 
resolution (number of gray levels) was set at 128 for bands 4, 

5, and 6 while band 7 was set at 64 so that the trace of the 
profile would not overlap the other bands. The acquisition of 
the data Wcis generated by driving the cursor along the transc:ct 
and recording the radiance values of every other pixel (a soft- 
ware design constraint) while in the profiling mode. 
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Essentially, the data thus acquired should correspond with the 
ridges and valleys along the transect, that is, each dip in 
radiance values should correspond to a transition from a south- 
east face to a northwest facing slope and, conversely, each rise 
in value should indicate a transition from a northwest face to 
a southeast facinr^ slope. To verify this observation, the 
radiance values were then correlated with the topography as shown 
on the uses maps. Figure 3 dramatically illustrates the large 
fluctuations of the radiance values with the terrain slopes and 
the changing forest reflectance along Transect B. 

The results of this preliminary step, besides demonstrating the 
above correlation, revealed that the inclination of the majority 
of the slopes was not great enough to produce a true shadow and 
that the observed phenomenon was indeed caused by a combination 
of terrain slope. Sun elevation, and bidirectional reflectance of 
surface features. To demonstrate how this phenomenon could 
affect a spectral classification, an unsupervised classification 
was made of the area containing Transect B. The classification 
consisted of eight nonoverlapping classes with the spectral bounds 
for each class selected at the option of the computex software. 

The results are shown in figure 5. The classes symbolized by I 
and " repre«’‘»nt the northwest facing slope while the southeast 
slopes are represented by the sumbols &, #, %, and $. Grasses 
and water were indicated by ( and ', respectively. The outlined 
areas in figure 5 represent the pseudoshadows caused by the Sun's 
low grazing angle. From this classification it is apparent that 
approximately 35 percent of the trace of the transect is affected 
by the terrain characteristics and low Sun angle. 

Based on the results of the above classification and a need to 
acquire more precise data, a series of slopes containing a near 
uniform canopy of mixed hardwood was selected along the trace of 
Transect B for mensuration of the radiance values. The Image 100 
system was used again to acquire the data v/hile employing the 


9 


4 ** Hlf l««l nf If I I • • •441 f Illllll • •• lll^•^ 444 ♦^^*>.• * I • • 1 • • •• l.ilie. 4 * M • t ? * I *M ! | 

1114411 .M • 44 '*-* 4 WjI I* I • 411 "4444 - 4 |M HI! .V .4 44 i ' 4 l ; • M * * I, H U * • i { ( * t $ t * LIL'. * M II 

44 444 * 41 t.l I • ^ 4441 1 I f I W * * f | | M ! M * I \ t * * * * * I I ! I I I * I I ' lt&h/. 44 * 4 44 * 44 ll**' 4 ||*' 4 l 4 t«.***iifKlt 4 t fl» 6 l 4 l ; I | I 
4II44* 4444444441 m*A II ttt .1 (*tl 4 “***Wf 4 M 4 * 4 lll 44 - 44 | '‘ 4 l ' % 4 M| I ^ I 44 1 4 **M 4 l* • 4 • 4 |•M•l^••••(f(<t••|lf^ 4 ll |•4 



^ * *M II 4 M*. . 1*41 M 4 I - •• 4 | II IvM4 ’ I l>.^* •• M 141 • 16444 * 

i ^ • I*.* II m I 4 * 44 l 4*4 > sl • I It I 4 |^.. 1^1 44 441 I M I i * 4 4 I ^E. * * t 44 4 I 


> I • 


III 


4 l 4 **l‘.M.\ll 4 llMA* 

( • V ( O I < • i * * M 

Kr*tf 

444461 

•i •tf_ 

111 I *i I 
1 1 f I I I I 
lilfIMIt'i < 

’ 4 |i l 

»tll 4 *ll** 44 f • 

<166114 M 4446 * i • 

,*t 4 **l'* 4 llll 6 f (*6 
/Ifi** ’ 4 I 66 IIIII 6 U* * 

•l 4 l 6 Ull 4 --r‘ 6 u** 

• M 4I6III4* 1 I -illti ( • • 

•. 6 \ 66 l 4 f !ll 4 || 4 l(((*MI|(, 

_ 16614 * I III ••4414 441 * <( * * *M 
II I 44 * • I 4441 lit M.M I 16 * I (i * P 

tl4*4*4l66N6i6tllll66r • I v« *11 , 

|•*I 6^|44 •• 4 *-*f-* 4 lH 1 ((* 6 III 44 

1416664 “ 

| 4 ||| 444||4 
*i** 44 lllJ 

••nil" 44 ini 66 M* ( r • I I - 4 l I I jilit 

4 ||| 4 *"" 44 * 44 |ll 6 *i'* 4 l!’’"M 444 l| I 4 i< 

• • * * * 6 1 1 1 1 1 4 II 1 6 ( ( u 1 1 • • • " • 4 • f U I 4 I \ *'*’•* ^TiTini 4 4 I 
1 61661 I 4 1666 16 I ( I I 6 * I 4 4 * 4111 4 . i f I 4 | i 4 M 4 ti n^ ^ I 444416 *. 
:* 4 *| ri* 4 l 4 ttS*l((* 4 l I *44444444 1 I If.F.KI^ ^ilM 6 .V.IM 

44 I* I M I * * • ' 441 * ( *641 44 446 * 6666 l i' ' / l I I I fl^^M 

il 4 **llMI •• 4 J‘.* * • * ♦• 4 * 4 * M • { U 4 M 4 I 6 * * * M.Tll 4 I I I i • 

441 1666 ^ 1 i. 111666 • * •••« 4 * 1 ! ' * I I m 46 * * 6 * ’ I . 6 1 1 
f *11 I 4 lii U 4 i 111166 * n ” 4 ‘ ** ' 116611 II I 446 * 6 < ' 16 * f 
44 tlimf 411 166^^66 i 44 * 4416 * * * ’ J I I ‘ * 4 ♦ 1 1 * ( • * * * 

4444416 I 


• 6 C 

• 6 (' 


C(((llt>» 6 ll *4166 
(C( 6 l. 6 l**f *646 


I I I f I I I I 4 lM l*. 4 * i 6 M I I t 6664 ' 4 i 4 | 44 ' * 4 4 * * ; i I I 4 * * i * I * 1 I 4 * i ' * * 1 1 1 6 ^ I 16 It 1 6 
*.*.' *^. 4 t%M 6 *&‘-IM V I I *41 II m II » 16 ..144111111 4 . 6 M 6 * * * 66 I 4 * 4 I 66 * 

44166*6661 III 14441 44 4 IMH 116 44*11 11111611*1 I * I I 4 11 * 44/6 I II * * * I ! *| 
,lt 14 l - »*4 444 * t :«-*** 44 * • 44 * 4 **“ 4 l 6 ^ I 4 1 1 4 4 4 ‘^ * 4 1 1 I 4 * * I * * 666 * * * * ’ " r 4 
»M|4*4Mltf4 *4H4444M'44H44N4|\*|| ./ 6664! 4l4****IM44|**l***"l^4il' 

t.i JIII4 fl“* ’ • a* “ • 444 • 44 • * 44| 4 ; '••4*i f 44 * 1 1 * 4 I " 4441 » ( • 4 Jl||*i*4 

-j 44|U •• *411 • • • 1 • • • *4! t • 4I..4? ! 1 *441 1 ! •?*4l4** * U • : r 


4 4 4 1^4 »'"444’44l|»4444l 
4|f J“LmIM4* 44C.664* 4441 
h' -I4l/^4* 1 *444 •|4*l 

'/ ! Jj I i->*i-4** * * 


Tf I oiitno'v 

' • • 4 * I * • 4 1 6 6 * U i < * * * 4 1 6 6 II ijf • " Y 

•II4*- M^*.^6* • I ‘iiM'-rti :if(T*:/inrif 

1 1 J *6666611 1 * O 64* - 16 1 * * r:;'4a||A 


4 I.. 4 ? ! 1 *441 1 ! '^'gli** * U • : 1 *4144 

t> l 4 ii^n 44 * 44 ll 4 U 4 | I 44 * 4 l 4 4 4 I ** 46 66 I 1 I 1 * H * 1 I * I I I * I * 1 4 U * * ! * * I I 4 4 
*4444 44*444411441 144414 * * 44416 * * 61 M 1 1 414 * 1 14611641 *IC ( * ** If Il 4 
44414 ^ 4 $ * * 444*44111444441111666 I 4 l 1 1 4 l 6 l 4 M 46661 * ill 1614 44 1444 * 
‘ 44444416661 4111416141 I 14 IIIM 14111111*44414 

1 **U 6 ^ 14 44 * * 441 * 1 14661* 1 1 464 * * 4 l 4 l *4444 

•• - ^** 1664 : If **| 4 ll* 4 *l- 4 |||** 4 * 4 | 

Transect B 46 * *41 1 •* 66 ^ 444 • 44 * 4 | 44 U 44 I 6 I 4 

^ 6 !. 6 n 4«|/.4 lk .*,4 ' • 46 * * 4 ! 14661*1 * 41666611144414 !* 

n*l 411^166 I » 0 '^J« 4 * ^^^ 6 l I 4 46 /. 6 I* 1 • 4 ^a 4 • I I 4 16 * *61661444 * *4 

Ji »ir . J«l • I I iijL «4 66 >l** 146661 I 1 *41411 114161 * I 44 MM 4 

‘/ 4 »I 6 .'^. f ^«4 .I 4 ** 4 f * * 64 ***l* * 6*6611 If* 444 4 
I.I 66 V 6 ,-. ^^ 1 * 166 I I 4 ** *411 44 41 1 4 * 4 * ! *** 4 | 4 * 

I • *|u?fjJu.vA I I I ,61 I 1444 * 44166661 * * I I * * ' f 4^61 

.|Tri 6 -.|>t|; r"^ ^ 16144 * 1 * * 4166 ’ Il 4 l * I I* ! * * 4 f M 44 

UlJ^l'MlIln « ^ 44 * 41166 - U**M 6 l*fl** I’l 4 

M • * • I H ‘ ! r yl 6 ll '-. 6 * I 4 ; 466 v *4 44 ^ fill 

■ « ifj 


iC/itIc*! I J 


1 14116* 


I W4I* 


• 466 I* 


i • I 


111 i<«4i4i66;/ 


4 4 41 1.66 4 
6 I 4 M. 6 * 66 I 4 ** " • mV'.I* 

A .661 II . 61 ! : 44 lJ 4 Kt;i 4 
M 4444 I * 1 '* 44 44 16 ' * 4 * 


tl 666 l **44 




f 1 1 • 


1 1 


61 11 ‘in 


I < I 




II 


.ill 


1 1 


1114441 
16661161 
**l**lll 

• * * 1 1 1 I 1 

44*1111! 

•Mil 466 
•I l**46.’ 
1*116611 
4 I 666 II I 
I 66 I 6 III 
14441116 

ir.666l 16 
11.661116 
1 16666 1 1 
k**6l4ll , 

•4*4411": V 

•*444*** 
••“4l***fl I 

r ••••••■Mi’ (( 


Legend ; 


I ■ 


Northwest facing 
slopes (hardwoods) 


I ! ' 


« 

$ 

% 

Cl 


* 6.66 I 414 166 


166* 1*1 
44444441 f 4* 4l 1 f I ’ * 4 
114**4444416*6 I*' 
• *44*-*4f Iif4’ 1* 

16/ 614411 11 1 4*444 
I ’ Ml * I * 6 I 4 I 4 A 4 4 I I 
6**4*1664 * 4 4 44 

5 16 i 44*..' in <14 4 4 44 4 

i » .1 I 
.11 
* * 141 
I-.66J, 


' II J H 4 I • 


1||4 Qi{ Ilf 4 
* 6*.66U4r 41 Uit 
« i ( * 6 I I 4 4 I 
*6* 1 . • |4l 
= • “ 1 * « 1 ’ 1 1 _ 

4IU6* * *66114441111 I 
'44446* * S ' 6666 l*l44HMrt4 . i 
*44441*^.44 U(ll4iM4U.tnu il 
" * 4 M 4 1 ; •. * 4 II 4 ■ U6 ' * * I 1 <1 d . I 4 

'*41666 • III! II 

•nil*! *4111* •"-4« *4llll *Mi 
I4II66I *41111 I*iln ‘l<.6** 1 6* 

<416114 ** 4 ; 44 . 44 I LUl 4 ^ 4 f ’ < ' I r I 

I * * 46 * * I * 4 4 4 * * ! m 4 ' * * 4 ; 4 * 1 I, I i I ‘ 6 * » * 6 
*466* *1444* *411 Jlt:f**<l’ ? '* 41 .1144.6. *6 

I666*6*444 4||IM|4' * * 4 i * * * tll>il4 tl 1 1 1 f 
66661*11444***4 144**** *•• .1 .|*M* *6111 

66144* * 4U4444I6I 4* 44414* * 4 Nltt4i( *6'.6* 

14*4*1 |i* 44 44414* *4444*416:11 16 (*l 446 

4*4"*4ll4! *441 *4* 1 *4444444 I Hf4l** * 1*4411 -‘-*^111116114 

6llll!**"4 44 44****4l444*l***ll4"44' *661 ' " 4l * I 1 I ' 16661 II II I I If* * 4lMi6l 16 1 1 44 * 46 ' * * * 6 * * ' 

• * *41 4* 1 114*4411* r 161 44 1 4** * * I ! * I I 4 I 1 1 II * * 4 4 I JJ I * 16.6 * * (611 II66II4* JJ---- 

• * 4* * * " I 1644* * • 44* 4161 1461414 41* 4M66M* * I 4 1 16 ' * J V I I * I * 4 4 * J * * «666.'. < *•< Ilflf6**l"*"4444 


Southeast facing 
slopes ( ha rdwoods ) 


Grasses 

Water 



< 1 1 6 * 


1*1 
466 
161 4 1666 


f4|l6ll4ll4* 


4* * 44* * * 461141 * r ( ' 4*41166111144 4 4*41164*4* *6ll4|44 II '! 416*6661 T J 1 J J U I 4 1 * * 416666* * * * 64f 1 1 1 1 1 I 444444 * J 1 1 


Scale 1:33 650 




Figure 5 .— An unsupervised classification utilizing eight 
spectral classes in the area of transect B. 



cursor in a tanjct acquisition mode in order to insure optimum 
positioning of the cursor. Duo to the size of the cursor 
(2 by 2 pixels) the apparent shadows had to bo larqe enough to 
accommodate the cursor uo that it could bo positioned totally 
within the targeted area. 

Figure 6 illustrates the 8i>octral res|x>nse from the selected 
targets representing the northwest and southeast slopes. As 
indicated by the resi)onso in Landsat bands 6 and 7, there is 
sufficient separation to form at least two major spectral groups 
or classes based on the slope and as*>oct angle. Of major 
significance is the wide variance in radiance within the two 
groups. This is especially pronounced in band 7 and to a much 
lesser extent in bands 5 and 6, while band 4 is the least 
affected. Those variances are not necessarily attributed to 
only the terrain, slope and aspect, and Sun angle, but may bo 
attributed to a certain amount of inconsistencies within the 
vegetation cover, such as stage of growth, tree height, density, 
and specie mix. 

7. COMPUTATIONS 

The second objective of this study focused on the development o*" 
a means for reducing the effects of toriain slope and Sun angle 
on Landsat radiance values. A simplified approach was desirable 
considering the time allocated for the task and the additional 
expense required of a more sophisticated effort which may or may 
not prove effective. 

The model, heroin discussed, reduces the variance in spectral 
response cause by the terrain sloptJ and Sun angle. This was 
achieved by normalizing all measured slopes along the transect 
to an incident angle of zero degrees, thus reducing the varia- 
tion caused by the slope. Figure 7 illustrates the pertinent 
illumination conditions. From this simplified geometry, it can 
bo determined that the radiance (C) , represented by the relative 
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Radiance value 





Figure 6.— Landsat spectral signatures representing mixed oak 
class and showing signature variation due to aspect angle. 
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zenith 



S ■ slope angle 
0 ■ pun's angle of incidence 
Cl ■ elevation angle of sun 
e ■ 90® - a - S 


Figure 7.- Illustration of simplified 
radiation geometry. 


radiance counts for each pixel is 


C ■ k cos e n 

where k and n are constants and 0 can be determined from the 
known values of Sun elevation (a) and terrain slope (s) by 

0 ■ 90* - a - 8 

where the southeast slopes are positive in sign and the north** 
west slopes negative. 

This procedure assumes that all radiance values are a linear 
function of cosine theta and that small inconsistencies may be 
attributed to inaccuracies in mensuration and/or the spectral 
inconsistencies of like vegetation. These small variations can 
be reduced to a minimum by linear regression and the constants 
k and n estimated (see appendix A) . 


Utilizing the constants k and n, as derived for the linear 
regression of Landsat band 7 and an estimated cosine 6, a set of 
corrected radiance values can be calculated which is independent 
of slope from the following relationship: 



Cij + k^ (1 - cos 6^) 


where i and j refer to the sensor bands and pixels, respectively. 
For the complete derivation of the above formulae refer to 
appendix B. 


Table 2 represents the preprocessing of the Landsat data by an 
empirical transformation having the form of the above equation 
and where the constants k and n (appendix A) were derived from 
a linear regression based on related spectral values and terrain 
slopes. 
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TABLE 1. - MSS REFLECTANCE IN RADIANCE VALUES SHOWING ONCORPECTED 
VERSUS CORRECTED DATA AND RESULTANT STANDARD DEVIATIONS FOR EACH 
LANDSAT MSS BAND. 
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8. DISCUSSION 


The claesif Ication of foreat types, as compared to croplands, 
is by no moans a simple process, especially since forests are 
indigenous to hilly and mountainous terrain where the slopes 
and aspect angle may cause phonological changes to the vegetation. 
This effect of steep terrain is compounded in the winter months 
when the Sun is low on the horizon. The resulting low light- 
level, coupled with sloping tftrrain has a tendency to increase 
the relative contrast between the illuminated slope and the 
opposite slope to the obvious detriment of the classification. 

In an attempt to solve this problem, initial thoughts may con- 
sider a simple clustering algorithm where one cluster may define 
the illuminated slope and the other cluster the opposite slope. 
However, all slopes exist, not just two, as there is a continuum 
of values. 

In the method presented here, an attempt is made to increase the 
accuracy of classifying forest areas by reducing the variation 
in spectral response due to terrain slope and Sun angle. In 
other words, by reducing the original Landsat signatures to a 
surface normal to the Sun's elevation, the large variation 
between the Sun-lit slope and the opposite slope to the Sun 
would be minimized. This effect is diagrammatical ly illustrated 
in figure 8, where the mean values of the corrected data for 
Transect B were plotted in addition to the minimum/maximum 
radiance values of the targets. The reduction in variance can 
easily be seen when compared to the uncorrected data in figure 9. 

It should be noted that the data for band 7 were used to calculate 
the new estimated cosine of theta; hence, all previous values 
have been reduced to one common value of a zero variance. Band 
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7 was selected to derive the values for theta because of the 
sensitive response to the angle of slope and the fact that it is 
least affected by haze. 

9. CONCLUSIONS AND RECOMMENDATIONS 

The work which has been accomplished thus far indicates that the 
model has definite promise. During the study several approxima* 
tions were deliberately made in order to provide a quick assess- 
inent as to the potential of the method. The authors are in 
full agreement that the method could no doubt be improved by 
considering all significant factors. However, this experiment 
provided preliminary indications of success by reducing signature 
variance which in effect is 4t contributary cause of misclassifi- 
cation. Additional data and testing are required to confirm 
the applicability of this method under various conditions. Based 
upon the foregoing study concerning the effects of terrain slope 
and low Sun angle on forest albedo, the following tenative 
conclusions can be made: 

e The Landsat spectral response varies considerably in 
areas of steep slopes and low Sun angles 

e The large variation of spectral response influences 
discrimination accuracies to the detriment of the 
overall classification 

e The spectral variation can be reduced by employing 
a simplistic illumination mode] to minimize the 
influence of topographic slope and low Sun angle 

It is recommended that future work should be organized as follows: 

a. Phase 1 — The plans for the immediate future should center 
on completing the remaining six transects so that the 
variable factors of vegetation types, aspect, and slope 
can be correlated to the mensurated radiance values. 
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This phase will also contribute to the development of 
Improved techniques for mensuration of radiance values 
and determination of terrain slope and aspect angle. 

b. Phase 2 — This phase will utilize the data and expertise 
acquired In Phase 1 to develop a cost effective working 
model applicable to forest Inventory needs. It Is 
recommended that Washington County, Missouri, be desig- 
nated us the initial test site due to availability of 
both Imagery and ground data for analysis and assessment 
of the working model. 
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APPENDIX A 


TRANSECT B TARGET GRAPHS 
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APPENDIX B 


A METHOD OF REDUCING CLASSIFICATION 
ERROR DUE TO SLOPE AND LOW SUN ANGLE 

The physical quantity that corresponds to the pixel radiance of a 
given band of Landsat data is the radiant power detected by the 
sensor. From theoretical considerations the pixel radiance C is 
linearly proportional to the cosine of the angle of incidence (6) 
of the Sun onto the area represented by the pixel. A simple 
formula for the radiance of a pixel area is 

C ■ k cos 0 + n (1) 

where k and n are constants. For the ith sensor band and the 
jth pixel 


Cj^j ■ k^ cos 0j + n^ (2) 

By employing the linear statistical model with error term i.e. 

^ij “ *^i 0j + n^ + (3) 

and an adequate sample of pixel measurements C^^ and 0^ (0^ is 
deduced from slope and aspect measurements) , the constrants k^^ 
and n^ can be estimated by linear regression. This enables 
subsequent estimation of the regression line radiance values 
c^j for any Sun angle 0^, i.e., 

Cij ■ k^ cos 0j + n^ (4) 

Assuming that svirface feature information is superimposed onto 
the 0 dependent regression line, the basis for the radiance modi- 
fication is to replace the pixel radiance value C^^^ with the 
expected value at normal incidence ^0 ® 0®), plus the deviation 
at 0j (see fig. B-1) to obtain: 
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Cj* - (kj + n^) ♦ 


where is the modified pixel radiance and 


^ij “ ^ij ■ ^’ij 


(6) 


Substituting for in equation (4) gives 

^ij “ ^ij *^i 

In order to circumvent the necessity for computing each 6^ from 
slope and aspect measurements, cos 6^ can be estimated from 
equation (4) using the radiance data from one of the bands 
(band 7 is chosen in this study) . An estimate of cos is 
given by 


(cos 0^) ■ (C^j - n^)/k^ (0) 

Substitution into equation (6) given 

^ij " ^ij + 

This procedure enables pixel radiance modification in three bands 
at the expense of one band. Three band classification can be 
performed on this modified data which is essentially independent 
of any slope effects other than true shadow. 



